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I. INTRODUCTION 



The behavior of atoms, molecules and ions placed in a strong magnetic field has attracted 
a significant attention during the last two decades (see, for example, the review papers [1- 
3]). It is motivated by both pure theoretical interest and by possible practical applications 
in astrophysics and solid state physics. Prom the point of theory, such studies would lead 
to a creation of a theory of atoms and molecules in a magnetic field similar to a standard 
atomic-molecular physics. In practice, even the basic elements of such a theory - a knowledge 
of the energy levels of the simplest Coulomb systems which can exist in a magnetic field - 
can be important for interpretation of the spectra of white dwarfs (where a surface magnetic 
field ranges in S fa 10^ — 10^ G) and neutron stars where a surface magnetic field varies in 
B 10^^ — 10^^ G, and can be even B 10^^ — 10^^ G for the case of magnetars. 

It was conjectured long ago [4, 5] that unusual chemical compounds can appear in a strong 
magnetic field. In particular, it was suggested by M Ruderman [5] and then developed by 
his followers (see [2] and references therein) that the presence of a strong magnetic field can 
lead to the formation of linear hydrogenic neutral molecules (linear chains) situated along 
magnetic lines. It was assumed that in the ground state all electrons are in the same spin 
state with all spins antiparallel to the magnetic field line. To avoid a contradiction with the 
Pauli principle it was further assumed that all electrons have different magnetic quantum 
numbers. It was considered as a characteristics of the ground state. It seems obviously 
correct in the case of atoms and atomic ions where the electrons are close to each other. 
However, it is not that obvious for the case of molecules where the electrons are situated 
in far distant places in space. All of them (or, at least, some of them) can be in the same 
quantum state, with the same spin projection and magnetic quantum number [6]. This 
situation was observed in H2 [7] and H3 [8], where in a domain of large magnetic fields 
the ground state was given by the state of the maximal total spin but with the electrons 
having the same zero magnetic quantum number (see a discussion below). In [5] qualitative 
arguments were presented that such chains can be of any length, thus, can contain arbitrary 
many protons. It seems that such a picture is oversimplified, it intrinsically assumes that 
the magnetic field is "infinitely" strong. Por instance, for any exotic chain (which does not 
exist in field-free case) there must be a certain threshold magnetic field since that it begins 
to exists. It can well happen that such a threshold magnetic field can be beyond of realistic 
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magnetic fields which occur in Nature. This phenomenon is absent in the quahtative theory 
[5]. Thus, some very general features of the Ruderman's picture only, like growth of the 
binding energies, shrinking of the size of the molecules with a magnetic field increase and 
maximal total electronic spin can hold for realistic high magnetic fields. 

It is well known that in absence of a magnetic field, in general, the hydrogenic linear 
chains (polymers) do not exist with the only exception of two shortest ones, H2 and H2 [20] . 
Therefore, for each other chain there must occur a threshold magnetic field from which the 
chain begins to exist if it is realized. It seems natural to assume that the threshold magnetic 
field grows with the length of the chain which is defined by a number of hescvy particles 
therein. At the moment, only those and H2 - the shortest chains - are studied in details, 
see e.g. [3] and [7] , respectively. The results are far more sophisticated than those predicted 
in a simple qualitative picture in [5]. For example, the H2-molecule does not exist at a large 
domain of strong magnetic fields. 

The aim of this article is to perform a detailed quantitative study of Hydrogen, He- 
lium and also mixed, Helium-Hydrogen linear chains with one or two electrons making an 
emphasis of the domain of magnetic fields 10^ < B < 10^a.u.(= 2.35 x 10^^ G). It is 
shown that in the one electron case depending on the magnetic field strength the hydro- 
genic systems H2 , Hg"*", Hi"*" and even Hg"*" can exist in linear geometry. It is also shown 
that, as the magnetic field is increased, the exotic helium-hydrogenic chains Heg"*", (HeH)^"*", 
(HHeH)^"^, (HeHHe)^"^ and Hel'^ begin to exist in linear geometry (see for a brief review 
[6]). For all magnetic fields the system Hg is stable when the system H3''" becomes stable at 
B > 10^^ G. A detailed review of the current status of some one-electron hydrogenic molec- 
ular systems, both traditional and exotic, that might exist in a magnetic field B > 10^ G 
can be found in [3]. For two-electron case depending on the magnetic field strength the 
hydrogenic chains H2, H3 , H4"'" and at most Hg"*" can exist in linear geometry, as well as the 
two-electron Helium chains He2''", Hcg"'" and He4''" , and the mixed Hydrogen-Helium chains 
(HeH)+, (H - He - H)2+,(He - H - He)=^+, (H - He - He - H)^+, (He - H - H - He)^+, 
(H-H-He-H-H)^+, (H - He - H - He - H)5+ and (He - He - H - He - He)^+. 
Since our study is limited to the question of the existence of a particular Coulomb sys- 
tem the main attention is paid to an exploration of the ground state. Overall study is made 
in framework of non-relativistic consideration by solving the Schroedinger equation. It is 
also assumed that the Born-Oppenheimer approximation of zero order holds, which implies 
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that the positions of positively-charged heavy particles are kept fixed (they are assumed 
to be infinitely- massive) . Relativistic corrections are always neglected assuming that the 
longitudinal motion of electrons is non- relativistic for magnetic field < 10^^ G while there 
are no relativistic corrections to the energies of transverse motion since the spectra of non- 
relativistic and relativistic harmonic oscillators coincide (we call it 'the Duncan argument', 
for a discussion see [9]). Some preliminary results were announced in [6]. 

Atomic units are used throughout {h—rrie—e—l), although energies are expressed in Ryd- 
bergs (Ry). The magnetic field B is given in a.u. with a conversion factor Bq — 2.35 x 10^ G. 

II. ONE-ELECTRON HYDROGENIC CHAINS 
A. Generalities 

Let us consider the electron and n infinitely-massive particles (protons) situated on a line 
which coincides to the magnetic line (sec Fig. 1). We call this system a linear finite chain 
of the size n. If for such a system a bound state can be found it implies the existence of the 
ion H^""^^''' in linear geometry. 



p p p p 



12 3 B n 

► 

FIG. 1: Hn" linear molecular ion in parallel configuration with a magnetic field B oriented 
along the z-axis. 

The Hamiltonian which describes this system when the magnetic field is oriented along 
the z direction, B = (0, 0, B) is [21] 

nn^{p + Ar-2j2-+ E ^ + 2B-S, (1) 

i=l,n i^j 

i,j=l,n 

(see Fig. 1 for the geometrical setting and notations), where — Zj — 1 in the case of 
protons, p = — iV is the momentum of the electron and S is the operator of the spin, is 



the distance from the electron to the ith proton and Rij is the distance between the ith and 
jth protons. .4 is a vector potential which corresponds to the constant uniform magnetic 
field B. It is chosen to be in the symmetric gauge, 

^=^(Bx r) = X, 0). (2) 

Finally, the Hamiltonian can be written as 

= (-v^ + - 2 f + ^ + + 24) , (3) 

^ ^ i=l * i^j 

where and Sz are the z-components of the total angular momentum and total spin 
operators, respectively, and p — -s/ + y^. Both and are integrals of motion. Thus, 
the operators and in (3) can be replaced by their eigenvalues m and respectively. 
Since we are interested by the ground state for which m — and = —1/2, the last term 
in (3) can be omitted and the reference point for energy becomes {—B). 

In the equilibrium configuration the problem is characterized by two integrals of motion: 
(i) angular momentum projection m on the magnetic field direction (2;-direction) and (ii) 
spatial parity p. The problem for parallel symmetric configuration is characterized by the 
2;-parity, Pz{z — > —z) with eigenvalues a — ±1. One can relate the magnetic quantum 
number m, spatial parity p and z-parity a, 

In the case m is even, both parities coincide, p — a. Thus, any eigenstate has two definite 
quantum numbers: the magnetic quantum number m and the parity p with respect f — > 
— r. Therefore the space of eigenstates is split into subspaces (sectors) each of them is 
characterized by definite m and a, or m and p. Notation for the states is based on the 
following convention: the first number corresponds to the number of excitation - "principal 
quantum number", e.g. the number 1 is assigned to the ground state, then a Greek letter 
a, IT, 5 corresponds to m = 0,-1,-2, respectively, with subscript g/u (gerade/ungerade) 
corresponding positive/negative eigenvalues of spacial parity operator P. 
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B. Method 



The variational method is used for a study of the Hamiltonian (3). Trial functions are 
chosen following the physics relevance arguments [10]. Their explicit expression is a linear 
superposition of K terms given by 

fc=l I ) k 

(see [3]), where and otk^i, h are linear and non- linear parameters, respectively. Inter- 
proton distances R are considered as variational parameters as well. Notation {} means 
the symmetrization of identical nuclei of the expression inside the brackets. Usually, to 
each term in (4) a certain physical meaning is given. For example, one term had all ctfe^j 
at i = 1, . . .n equal being an analogue of the Heitler-London wavefunction for the H^-ion 
- describing the coherent interaction of the electron with all protons. For another term all 
ctfc^i, except for one, vanish being an analogue of the Hund-Mulliken wavefunction - describ- 
ing the incoherent interaction of the electron with all protons. All other terms are different 
non-linear superposition of these two being an analogue of Guillemin-Zener wavefunction 
for the Hj-ion. We call a term for which all a^^i are different and unconstrained, the general 
term. Needless to mention that in each particular term in (4) the parameters are chosen in 
such a way to assure normalizability of this term as the overall function. 

Calculations are performed using the minimization package MINUIT from CERN-LIB. 
Two-dimensional integration is carried out using a dynamical partitioning procedure: a 
domain of integration is manually divided into subdomains following an integrand profile 
with a localization of domains of large gradients of the integrand. Each subdomain is 
integrated (for details, see, e.g., [3]). Numerical integration of subdomains is done with a 
relative accuracy of ~ 10~^ — 10~^° by use of the adaptive DOIFCF routine from NAG-LIB. 

n = 1. This case was considered for the sake of completeness. It is known that the hydrogen 
atom exists for any magnetic field strength. It is the least bound system among one- 
electron systems. The results for H-atom B — 10^, 10^ a.u. are calculated with 
a ten-parametric variational trial function which is a modification of the function 
introduced in [10, 11]. It will be described elsewhere. 

n = 2. The results for H^-ion are found with 3-term trial function (4) which depends on the 10 
free parameters including the interproton distance it is a linear superposition of the 
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Heitler-London, Hund-MuUiken and Guillemin-Zener (general term) wavef unctions. 
For B < 10^ a.u. results are from [3]. 

n = 3. The results for H3'''-ion are found with a 3-term trial function (4) which depends on 22 
free parameters including two interproton distances i?'s, it is a linear superposition of 
the Heitler-London, Hund-MuUiken and a type of the Guillemin-Zener (general term) 
wavefunctions. For B < 10^ a.u. results are from [3]. 

n = 4. Results for H4"^-ion are found with 1-term trial function (4) which depends on the 7 
free parameters including three interproton distances i?'s two of them are assumed to 
be equal (symmetric configuration). For B < 10^ a.u. the results obtained with 3- 
and 7-term trial function (4) can be found in [3]. They lead to slightly better binding 
energies but do not change the qualitative picture. 

n = 5. It is the first study of this system. The results for the H5"'"-ion are obtained using 
a 2-term trial function (4) which depends on the 15 free parameters including four 
interproton distances i?'s, two pairs of them are assumed to be equal (symmetric 
configuration). In fact, it implies that a linear superposition of two general terms is 
taken. It is worth noting that already 1-term trial function at B — 10^ a.u. gives a 
clear indication to the existence of the Hg'^-ion with binding energy Ej, — 206.11 Ry 
and equilibrium distances Ri = 0.053 a.u., R2 = 0.032 a.u. The smallest magnetic field 
for which a minimum of the total energy surface in i?'s was observed is 5 x 10® a.u. The 
H5"'"-ion for these magnetic fields looks like H3"'"-ion bound with a far-distant proton 
from each side. 

n = 6. No indication to the existence of the Hg^-ion in the domain B < 10^ a.u. is found. 
C. Results 

The results of the calculations are presented in Tables I-II. Two traditional for field-free 
case systems H and H2 exist for all studied magnetic fields B < 10^ a.u. The first exotic 
molecular system Hg^ appears at B ^ 10^ a.u. and exists for larger magnetic fields. Another 
exotic molecular system appears at -B ~ 10"^ a.u. and the last exotic molecular system 
appears at -B ~ 5 x 10® a.u. No other one-electron molecular hydrogenic systems are 
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seen for B < lO^a.u. For n > 1 the optimal geometry of any molecular system is linear and 
aligned along magnetic field. Thus, such a system forms a finite chain. It is checked that 
the configuration is stable with respect to small deviations from linearity. All studied finite 
chains arc characterized by two features: with a magnetic field growth (i) their total energies 
increase and (ii) their lengths decrease - each system becomes more bound and compact. 

For all studied magnetic fields the systems H and are stable: the H-atom has no decay 
channels, where the total energy of the Hj-ion is always less than the total energy of the 
H-atom. Furthermore, for B < 1.5 x lO^a.u. the H^-ion has the smaller total energy then 
Hg^'^-ion when exists: these two finite chains are the only ones which exists in this domain. 
The Hg^-ion never dissociates to H + 2p but it always dissociates to + p. For higher 
magnetic fields B > 1.5 x 10^ a. u. the Hg^-ion becomes stable as well. It is characterized 
by the smallest total energy for these magnetic fields. Another exotic molecular system 
114"*" never dissociates to H + 3p, but it dissociates to + p for 10^ < B < 10^ a.u. For 
magnetic fields B > 10^ a.u. the total energy of 114"^ is smaller than and the latter 
dissociation channel does not exist. For all studied magnetic fields B < 10^ a.u. the system 
114"*" can dissociate to Hg"*", although the energy difference between such systems decreases 
gradually as the magnetic field increases. A smooth extrapolation indicates that at the 
magnetic B 2 x 10^ a.u. there is a crossing for which the total energies of Hg"*" and 114"*" 
become equal. The system Hg"*" can dissociate to all finite chains except for single proton 
one: H-atom. Summarizing, one can state that there are two one-electron finite hydrogenic 
chains characterized by lowest total energy for different magnetic fields: it is the H^-system 
at < 5 < 1.5 X 10^ a.u. and the Yl\^-ion at 1.5 x 10^ < 5 < 10^ a.u. 

III. TWO-ELECTRON HYDROGENIC CHAINS 
A. Generalities 

Let us consider a system of two electrons and n infinitely-massive protons situated on 
a line which coincides to the magnetic fine (see Fig. 2). It is called 2e-linear finite chain 
of the size n. If for such a system a bound state can be found it implies the existence of 
the ion H^T"^^"^ in hnear geometry. Sometimes, we say that above system is "in the parallel 
configuration". Also, it implies that the corresponding finite chain exists. It can be stable 
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B(a.u.) 

System 


1 10 10^ 10^ 10^ 10^ 


H 


1.0 1.662 3.495 7.564 27.10 73.96 108.86 




1.2053 1.9499 4.3498 10.291 45.799 139.91 217.75 




- - - 8.639 45.408 160.17 263.80 




- - - - 34.922 142.75 251.71 




------ 206.15 



TABLE I: Binding energies (in Ry) for the ground state lag of the one-electron hydrogenic hnear 
systems (finite chains) in a magnetic field. Binding energies for the ground state Isq of the H-atom 
at < 5 < 10^ a.u. from [11]. 



B(a.u.) 

System 





1 


10 


102 


10"^ 


10^ 


10^ 


H+ (R) 


1.997 


1.752 


0.957 


0.448 


0.118 


0.045 


0.032 


Rl+ {R,R) 








0.579, - 


0.130, - 


0.044, - 


0.029, - 


hI"*" {Ri,R2,Ri) 










0.214, 0.138, - 


0.056, 0.044 , - 


0.034 , 0.028 , - 


{Ri,R2,R2,Ri) 














0.053, 0.032, - , - 



TABLE IL Interproton equilibrium distances (in a.u.) for the ground state lag of the one-electron 
hydrogenic linear systems (finite chains) in a strong magnetic field. All configurations have center 
of symmetry, symmetric interproton distances are not displayed. 

or metastable. 

The Hamiltonian which describes the system of two electrons and n protons when the 
magnetic field is oriented along the z direction, B = (0, 0, B) is [19] 
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B 
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FIG. 2: H^" ^■'^ linear molecular ion in parallel configuration with a magnetic field B oriented 
along the z-axis. 



W„ = E(P' + -4<)'-2Er + f + E ¥^ + 2B.S, (5) 

i=l,n i,j=l,n 

(see Fig. 2 for the geometrical setting and notations), where Zi — Zj — 1 in the case of 
protons, Pi = —iVe is the momentum of the £th electron, ru is the distance from the ^th 
electron to the ith proton and Rij is the distance between ith and jth proton, = \ fi —f2\ 
is the interelectron distance, where fl (r^) is the position from the center of the chain (mid- 
point with respect to the end-situated protons) of the first (second) electron and S = S1-I-S2 
is the operator of the total spin. Ai is a vector potential which corresponds to the constant 
uniform magnetic field B written in the symmetric gauge (2). Finally, the Hamiltonian can 
be written as 

i=l,n i,j=l,n 

where — L^^+Lz^ and Sz — Sz^+Sz^ are the z-components of the total angular momentum 
and total spin, respectively, and = \/ x\ + yj. All performed calculations released a 
symmetry property of a chain: in the optimal geometry a chain has a center of symmetry. 
Hence, for any proton there is a partner situated symmetrically with respect to this center. 
We consider that property as intrinsic of any chain. 

The problem under study is characterized by three conserved quantities: (i) the operator 
of the 2;-component of the total angular momentum (projection of the angular momentum on 
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the magnetic field direction) giving rise to the magnetic quantum number m, (ii) the spatial 
parity operator P{fi — > —fi,f2 —^2) which has eigenvalues p — ±l(gerade/ungerade) 
(iii) the operator of the 2;-component of the total spin (projection of the total spin on the 
magnetic field direction) giving rise to the total spin projection m^. Hence, any eigenstate 
has three explicit quantum numbers assigned: the magnetic quantum number m, the total 
spin projection and the parity p. For the case of two electrons the total spin projection 
nis takes values 0, ±1. 

As a magnetic field increases a contribution from the Zeeman term (interaction of spin 
with magnetic field, B • S) becomes more and more important. It seems natural to assume 
that for small magnetic fields a spin-singlet state is the state of lowest total energy, while 
for larger magnetic fields it should be a spin-triplet state with = — 1, where the electron 
spins are antiparallel to the magnetic field direction B. The total space of eigenstates is 
split into subspaces (sectors), each of them is characterized by definite values of m, p and 
777.5. It is worth noting that the Hamiltonian Hn is invariant with respect to refiections P^: 
2^1 —> —^1 and Z2 — > —Z2, with eigenvalues crjv = ±1, for a symmetric chain. 

In order to classify eigenstates we follow the convention widely accepted in molecular 
physics using the quantum numbers m,p and the total spin S without indication to the 
value of nis- Eventually, the notation is '^^~^^Mp, where 25-1-1 is the spin multiplicity which 
is equal to 1 for spin-singlet state {S — 0) and 3 for spin-triplet {S — 1), as for the label 
M we use Greek letters E,n, A that mark the states with \m\ — 0,1,2,..., respectively, 
but implying that 777 takes negative values and the subscript p (the spatial parity quantum 
number) takes gerade/ungeradc((y'/M) labels describing positive p = +1 and negative p = —1 
parity, respectively. There exists a relation between the quantum numbers corresponding to 
the ^-parity and the spatial parity: 



Present consideration is limited to the states with magnetic quantum numbers m — 0,-1 
because the total energy of the lowest energy state in any sector with positive 777 > is 
always larger than one with 771 < 0. A study of states with different m is necessary to 
identify the state of lowest total energy. At large magnetic fields for all studied two-electron 
chains this state was characterized by 777 = —1 in agreement with Ruderman's hypothesis. 
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B. Method 



As a method to explore the problem we use the variational procedure. The recipe of 
choice of trial functions is based on physical arguments [10]. As a result the trial function 
for the lowest energy state with magnetic quantum number m is chosen in the form 

fe=l I ) k 

where Ue — ±1 stands for spin singlet (+) and triplet states (— ), while {} means the 
symmetrization of identical nuclei of the expression inside the brackets. The P12 is the per- 
mutation operator for the electrons, (1 -^-^ 2). The a^^ij, Pk,i-2 and 7^ as well as interproton 
distances Rij = Rji are variational parameters. For each term with fixed k their total num- 
ber is 2fi + 4 including the linear parameter A^. In addition, we have n — 1 interproton 
distances. It is worth emphasizing that in the trial function (7) the interelectron interaction 
is included explicitly in the exponential form e'^^^'^. 

Calculations are performed using the minimization package MINUIT from CERN-LIB. 
Multidimensional integration is carried out using a dynamical partitioning procedure: a 
domain of integration is manually divided into subdomains following an integrand profile 
with a localization of domains of large gradients of the integrand. Each subdomain is 
integrated separately using parallelization procedure (for details, see, e.g., [3]). Numerical 
integration of subdomains is done with a relative accuracy of ~ 10~^ — 10~^ by use of the 
adaptive DOIFCF routine from NAG- LIB. A process of minimization for each given magnetic 
field and for any particular state was quite time-consuming due to a complicated profile of 
the total energy surface in the parameter space but when a minimum is found it takes several 
seconds of CPU time to compute a variational energy. 

n = 1. This case corresponds to the negative hydrogen ion H~ and is mentioned for the sake of 
completeness. It is known that the negative hydrogen ion H~ exists for any magnetic 
field strength [12]. At zero and small magnetic fields i? < 5 x 10^^ a.u. the spin- 
singlet state ^0 is the ground state, li B > 5 x 10~^ a.u. the spin-triplet state ^(—1) 
which does not exist in the absence of a magnetic field becomes bound and the ground 
state. Although this result is checked quantitatively for magnetic fields up to 4000 a.u. 
[12, 13] it is quite hkely that it holds for higher magnetic fields. It is the least bound 
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5(a.u.) 
System 


10^ 10^ 10^ 


4.414 X 10^=^ G 


10^ 


10^ 




8.35 16.95 30.1 


35.4 


82.5 


121.4 


H2 


16.473^ 35.632 71.42 


85.00 


219.9 


330.3 




18.915 44.538 95.21 


115.19 


324.2 


529.8 




(17.601) (43.917) 99.80 


122.34 


367.7 


636.0 




91.70 


114.34 


383.2 


687.7 



TABLE III: Double ionization energies Ej in Ry {Et = —Ej) for the ground state of the 
two-electron hydrogenic systems (finite chains) in a strong magnetic field. ^ from [7]. Energy in 
brackets means that the state is bound but the ground state corresponds to an unbound state. 
The magnetic field Bschwinger = 4.414 x 10^^ G = 1.878 x 10^ a.u. corresponds to the so called 
non-relativistic threshold for which the electron cyclotron energy equals the electron rest mass. 

system among two-electron systems made from protons. However, the H~-ion is stable 
for studied magnetic fields: the dissociation H~ — > H -|- e is prohibited. 

n = 2. In a domain of non ultra- high magnetic fields the H2-molecule was studied in details 
in [7]. It was shown that the lowest total energy state depends on the magnetic field 
strength. It evolves from the spin-singlet ^T,g state at < S < 0.18 a.u. to a repulsive 
spin-triplet state (unbound state) for 0.18 a.u. ^ B < 12.3 a.u. and, finally, to 
a strongly bound spin-triplet state. Hence, there exists quite large domain of 
magnetic fields where the H2-molecule is unbound being represented by two hydrogen 
atoms in the same electron spin state but situated at infinite distance from each other. 
The optimal geometry of the H2-molecule (when exists) corresponds always to the 
elongation along a magnetic line for the ^E^, state thus forming a finite chain. It is 
assumed that the chain in the ^H„ state is stable towards the deviation from linearity. 
This assumption seems well justified from physics point of view at large magnetic 
fields: any deviation from linearity leads to a sharp increase in total energy due to 
non-vanishing rotational energy. This chain is stable (when exists) for all studied 
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B{a.u.) 

System 


10^ 


10^ 


10^ 


4.414 X 10^3 G 


10^ 


10^ 


H2 (R) 
H3 {R, R) 

{Ri,R2,Ri) 

Hf"*" i?2, -Ri) 


0.38^ 
0.395, - 
(0.51, 0.38, -) 


0.19 
0.183, - 
(0.215, 0.174, -) 


0.102 
0.093, - 
0.103, 0.086, - 
0.184,0.134 , - , - 


0.087 
0.078, - 
0.092, 0.075, - 
0.160,0.110 , - , - 


0.038 
0.030, - 
0.030, 0.018 , - 
0.035,0.025,-,- 


0.034 
0.023, - 
0.020, 0.013, - 
0.023 , 0.018 , - , - 



TABLE IV: Interproton equilibrium distances (in a.u.) for the ground state ^n^j of the two-electron hydrogenic linear systems (finite chains) 
in a strong magnetic field. All configurations have center of symmetry and symmetric interproton distances are not displayed. * from [7]. 
Distances in brackets mean that the state is bound but the ground state corresponds to an unbound state. 



magnetic fields. However, this chain always has the total energy higher than the Hg 
chain (see below) and thus less preferable energetically. Calculations for the state 
of H2 using a single function of the form (7) for which all a parameters are different 
(general term) are presented in Tables III-IV. 

3. In [8] it is shown that the H3 molecular ion exists in a magnetic field as a bound state. 
For B > 0.2 a.u. the ground state geometry is realized in the linear, parallel to the 
magnetic field line configuration. Thus, the three-proton finite chain occurs. In the 
domain 0.2 ^ B < 20 a.u. the ground state is realized by state and it is weakly 
bound. However, at > 20 a.u. the ground state ^H^ state is strongly bound and the 
chain is stable. 

4. In field-free case the system (4p2e) does not display any binding. However, for mag- 
netic fields B > 2000 a.u. it becomes bound in the linear configuration aligned along 
the magnetic fine with the ^H„ state as the ground state. Hence, the molecular ion 
H4''' begins to exist. Its total energy is lower systematically than the total energy of 
HJ. Hence, the molecular ion H^+ is stable. With an increase of the magnetic field 
strength, the total energy at the equilibrium position decreases, the system becomes 
more bound (in this case, the double ionization energy is Ej — — -E^T; it increases 
with B) and more compact (the internuclear equilibrium distance decreases with B). 
Eventually, we state that the finite chain R^'^ is always stable. For magnetic fields 
1 < B < 2000 a.u. the state ^H^ is bound but the ground state corresponds to an un- 
bound system in the repulsive state: It consists of two H2 ions at infinite distance 
from each other. 

5. In field-free case the system (5p2e) does not display any binding. However, for mag- 
netic fields B > 5000 a.u. it becomes bound in the linear configuration aligned along 
the magnetic line with the state as the ground state; hence, the molecular ion Hg^ 
begins to exist. For 5000 ^ B < 10^ a.u. the Hg^ molecular ion decays to H4+ -|- p. 
At magnetic fields B > 10^ a.u. the molecular ion H5+ becomes stable. 

6. It is not found an indication to the bound state of the (6p2e)-system even for the 
highest magnetic field studied. 
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C. Results 



The results of the calculations are presented in Tables III-IV. Three traditional for 
the field- free case systems H~, H2 and continue to exist at magnetic fields lO^a.u. < 
B < lO^a.u. The first exotic molecular system H4''" appears at ~ 2 x lO^a.u. in linear 
configuration and exists for all larger magnetic fields. Another exotic molecular system Hg"*" 
appears at slightly larger magnetic field ~ 5 x 10^ a.u. No other two-clcctron molecular 
hydrogenic system is seen in the domain B < 10^ a.u. At large magnetic fields the ground 
state of all studied systems is the spin-triplet state with spin projection mg — —1 and total 
magnetic quantum number m = — 1 in agreement with Ruderman hypothesis. For n > 1 
the optimal geometry of the molecular system is linear, and the system is aligned along 
magnetic field. Thus, each molecular system forms a finite chain. It is checked that such a 
linear configuration is stable with respect to small vibrations and its vibrational energies can 
be calculated. However, we were not able to check stability of the configuration with respect 
to small deviations from linearity and to calculate the rotational energies. All studied finite 
chains are characterized by two features: with a magnetic field growth (i) their binding 
energies increase and (ii) their longitudinal lengths decrease - each system becomes more 
bound and compact. For all studied magnetic fields B > 10^ a.u. the systems H~ and 
H2 arc stable. They are characterized by much smaller binding energies in comparison with 
other systems. Thus, their significance for a thermodynamics at a fixed magnetic field seems 
limited. 

It is worth emphasizing that among two-electron hydrogenic finite chains the system H3 
has the lowest total energy in the domain 10^ < B <2 x 10^ a.u., at larger magnetic fields 
2 X 10^ ^ B < 10^ a.u. the finite chain 114"*" gets the lowest total energy and, eventually, at 
B > 10® a.u. the molecular ion H5''' (the longest hydrogenic chain) is characterized by the 
lowest total energy. Interestingly, in the domain 10® ^ B < 10^ a.u. all two-electron finite 
Hydrogen chains are stable. 
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IV. ONE-ELECTRON HELIUM AND HELIUM-HYDROGEN CHAINS 



A. Generalities 

Let us consider now molecular systems composed of one electron and a finite number n 
of infinitely- massive protons and/or a-particles as charged centers, situated on a line which 
coincides to the direction of an homogeneous magnetic field. The geometrical arrangement 
is similar to that depicted in Figure 1, except for the fact that charged centers can be 
either protons or a particles. If found, bound states of such systems are called one-electron 
helium or helium-hydrogen chains. In the present review only one-electron helium or helium- 
hydrogen chains with n — 1,2,3 were included. 

Following similar considerations as for the case of hydrogenic chains (see section II), 
the Hamiltonian which describes the one-electron helium (helium-hydrogen) chains in a 
magnetic field oriented along the z direction, B = (0, 0, B) is given by the Hamiltonian (3) 
with Zi, Zj = 1 or 2, depending on each particular system. Since we are interested by the 
ground state for which m = and = —1/2, the last term in (3) can be omitted and the 
reference point for energy becomes equal to {—B). 

B. Method 

The variational method is used for a study of the helium (helium-hydrogen) chains de- 
scribed by the Hamiltonian (3). Trial functions arc chosen following physics relevance argu- 
ments [10]. Their explicit expressions are linear superpositions of K terms given by functions 
of the class (4), where and a^^i, (3k are linear and non- linear parameters, respectively. In- 
ternuclear distances R are considered as variational parameters as well. In this case the 
notation {} in (4) means the symmetrization of the expression inside the brackets with 
respect to the permutations of the identical charged centers (for example for the system 
(HHeH)^+ it means permutation with respect to the external protons. As for the case of 
hydrogenic chains, each term in (4) has a certain physical meaning (see section II). In the 
following we describe the different chains studied. 

n = 1. (ae). This case corresponds to the simplest one electron helium system. It is known 
that the positive atomic ion of helium exists for any magnetic field strength. Fur- 
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thermore, it is the only one electron helium system which exists for magnetic fields 
of strength B < 10 a.u. The results presented below for the ground state Isq of the 
He"*" atomic ion (see Table V) were obtained with a seven-parametric variational trial 
function introduced in [11] for a study of the H-atom. 

n = 2. (i) (aae). Accurate variational calculations in equilibrium configuration (parallel to 
the magnetic field) for the ground state lag of the system Re^'^ were carried out 
in details in [3, 17] for the range of magnetic fields 10^ a.u. < -B < -Bschwinger- A 
3-term trial function of the form (4) which depends on ten free parameters includ- 
ing the internuclear distance R is used in the calculations. It is the same linear 
superposition of the Heitler-London, Hund-MuUiken and Guillemin-Zener wave- 
functions which was used to study the molecular ion (see section II above) . It 
is found that for magnetic fields 10^ B < 10^ a.u. the system He2''' is unstable 
towards the decay to He"^ + a. Nonetheless, at B > 10'^ a.u. this compound 
becomes the system with the lowest total energy among the one electron helium 
(helium-hydrogen) chains. In [17] lowest vibrational and rotational energies for 
this system were also calculated. 

(ii) (ape). The first indication about the existence of the hybrid system (HeH)^'*', for 
magnetic fields B > 10^ a.u., was established in [3, 17], where accurate variational 
calculations for the ground state la of the system (HeH)^"*" were carried out. 
Variational calculations are done with a 3-term trial function of the type (4) . In 
[3, 17] it was also demonstrated that the equilibrium configuration corresponds 
to the situation when the molecular axis (the line connecting the proton and the 
a particle) is parallel to the magnetic field. For the narrow range of magnetic 
fields 10"^ a.u. ^ B < Bschwinger the system (HcH)^^ is found to be a long-living 
metastable state decaying to He'*' -|- p. For magnetic fields larger than Bschwinger 
the system becomes stable towards the decay to He"*" -|- p. 

n = 3. (i) (aaae). It seems it is for the first time wc see an indication to the possible 
existence of the exotic molecular ion Hcg'*' for magnetic fields B > 10^ a.u. For 
this system a 3-term trial function of the form (4) is used for its variational study. 
It depends on 22 free parameters including two internuclear distances Ri^2- This 
function is the same linear superposition of the Heitler-London, Hund-MuUiken 
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and a type of the Guillemin-Zener wavefunctions which was used to study the Hg"*" 
molecular ion (see section II above). It is found that the system {aaae) begins 
to exist as a bound state (i.e. displays a minimum in the corresponding potential 

energy surface for finite internuclear distances) at magnetic fields B > 10^ a.u. in 
the linear symmetric configuration (for which the two internuclear distances are 
equal, Ri — R2) parallel to the magnetic field direction. Ground state is lUg. 

(ii) (pape). First indications on the existence of the exotic trilinear molecular ion 
(H — He — H)^"*" for magnetic fields B > -Bschwinger, were given in [3, 17]. For 
this system a 3-term trial function of the form (4) which depends on 14 free 
parameters including two Ri^2 is used in the variational calculations. The results 
clearly show the appearance of a minimum in the potential energy surface of the 
(appe) system for the symmetric configuration of the charged centers (p — a — p) 
with R1 — R2. Ground state is the type lUg. It was not seen an indication to the 
existence of non-symmetric configuration (a — p — p). 

(iii) {apae). First indications on the existence of the exotic trilinear molecular ion 
(He — H — He)^"'' for magnetic fields B > -Bschwinger were given in [3, 17]. For this 
system a 3-term trial function (4) which depends on 14 free parameters including 
two internuclear distances i?i,i?2 is used in the variational calculations. The 

results show the appearance of a minimum in the potential energy surface of the 
(aape) system for the symmetric configuration of the charged centers (a — p — a) 
with R1 — R2. Ground state is lag. It was not seen an indication to the existence 
of non-symmetric configuration {a — a — p). 

n = 4. No binding is detected for systems (aaaae), (appae), (paape) even for the highest 
studied magnetic field ~ lO^a.u. 

C. Results 

The results of the ground state calculations are presented in Tables V-VI. The positive 
atomic ion of helium He^ is the only system which exists for all studied magnetic fields 
< -B < 10^ a.u. At i? ~ 10^ a.u. the first exotic molecular system Hc2^ appears being 
unstable towards decay to He'''-|-Q! in the range of magnetic fields 10^ a.u. < B < 2 x 10^ a.u. 
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For larger magnetic fields B > 2 x lO^a.u. the system He2"^ becomes the most bound 
one-electron system among the systems made out from protons and/or a-particles and it 
is stable. Two exotic molecular systems begin to exist at about the same magnetic field 
B ~ 10"^ a. u. Namely, the hybrid molecular ion (HeH)^^, followed by the trilinear symmetric 
molecular system (H — He — H)^^, being unstable towards decay to He^ + p and He^ + 2p, 
respectively. Remarkably, the system (HeH)^"'' rapidly becomes stable for magnetic fields 
B >2 X 10^ a.u. The system (H — He — H)^"^ becomes more bound than He"*" for magnetic 
fields B > 5 X 10^ a.u. but remains unstable towards decay to (HeH)^"^ in the range 
of magnetic fields 5 x 10^ < -B < 10^ a.u. It never dissociates to H^ + a. Another exotic 
symmetric molecular system (He — H — He)^"*" appears at 5 ~ -Bschwingcr, being unstable 
towards decay to (HeH)^'*' + a for magnetic fields -Bschwingcr ^i?<6.5xl0^ a.u., as well as 
towards decay to He2^ + p for all magnetic fields studied. It is worth noting that, in spite 
of the greater Coulomb repulsion, the system (He — H — He)^"*" becomes more bound than 
(H — He — H)^"^ for magnetic fields S > 1.8 x 10^ a.u. The last exotic molecular system Hcg"*" 
appears at B > 10^ a.u. This system is unstable with respect to decay into He2"^-|-a. Present 
level of available computational resources does allow to draw a reliable conclusion about this 
molecular system at larger magnetic fields. No more one-electron helium-hydrogenic system 
is seen for the range of magnetic fields studied B < 10^ a.u. 

Concrete variational calculations for the chains Hel'^ and (HeH)^"*" demonstrate that the 
optimal geometry of the molecular systems is linear and aligned along magnetic field, being 
stable with respect to small deviations from linearity. This is understood with simple argu- 
ments since any slight deviation from the magnetic field direction leads to a large increase 
in the rotational energy. So, it is natural to assume that all other studied linear chains are 
also stable with respect to small deviations from linearity. 

All studied finite chains are characterized by two features: with a magnetic field growth 
(i) their total energies increase and (ii) their equilibrium size decreases - each system becomes 
more bound and compact. 

Summarizing, one can state that among the one-electron helium-hydrogen chains there 
arc two helium systems characterized by the lowest total energy for different magnetic fields: 
it is the He+ ion at < B < 2 x 10^ a.u. and the Hco^-chain at 2 x 10^ < 5 < 10^ a.u. 
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B (a.u.) 
System 


1 


10 


102 


10^ 


4.414 X 10^3 G 


10^ 


10^ 


He+ 


4.8820 


8.7801 


19.109 


78.426 


92.528 


226.66 


345.17 








16.516 


86.233 


105.121 


305.11 


507.31 


Hel+ 












227.83 


417.15 


/TT TT\2-|- 

(HeH)^+ 








77.303 


92.858 


251.32 


402.10 


(HHeH)^+ 








64.747 


79.69 


233.71 


392.47 


(HeHHe)^+ 










70.76 


230.38 


408.58 



TABLE V: Binding energies (in Ry) for the ground state lag of the one-electron hchum and 
hehum-hydrogenic hnear systems (finite chains) in a magnetic field (the ground state for (HeH)^"'" 
is la). For He|+: Eb = 86.76 Ry, Req = 0.202 a.u. aX B = 10^^ G, while for Bschwinger there is no 
minimum. 

V. TWO-ELECTRON HELIUM AND HELIUM-HYDROGEN CHAINS 
A. Generalities 

Let us consider systems of two electrons and n infinitely-massive a-particles situated on 
a line which coincides to the magnetic line. If a bound state is found the system is called 
2e-linear Helium chain of the length n indicating the existence of the ion 
geometry. 

The Hamiltonian which describes systems of two electrons and a number of a particles 
with a magnetic field oriented along the z direction, B = (0, 0, B) is given by the Hamiltonian 
(6) with Zi = Zj = 2. 

All performed calculations show that in the optimal geometry the chain possesses a 
symmetry property similar to two-electron hydrogenic chains: for any a-particle there is a 
partner situated symmetrically with respect to the center of the chain. 

Another type of systems we study are mixed ones: out of n heavy centers some of them 
have the charge two (a-particles) and some have the charge one (protons). If a bound state 
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B (a.u.) 

System 


102 


10^ 


4.414 X 10^3 G 


10^ 


10^ 


He3+ (R) 


0.779 


0.150 


0.126 


0.049 


0.032 


He|+ {R, R) 








0.070, - 


0.041, - 


/XT TT\ V— I— / T-l\ 

(11611)"^+ (i?) 




0.142 


0.119 


0.048 


0.032 


(HHeH)^+ {R,R) 




0.227, - 


0.184, - 


0.058, - 


0.035, - 


(HeHHe)^+ {R, R) 






0.170, - 


0.051, - 


0.031, - 



TABLE VI: Internuclear equilibrium distances (in a.u.) for the ground state lag of the one- 
electron helium and helium-hydrogenic linear systems (finite chains) in a strong magnetic field 
(the ground state for (HeH)^''' is la). For all configurations which have center of symmetry, 
symmetric internuclear distances are not displayed. 

is found the system is called 2e-linear Helium-Hydrogen chain of the length n. 
B. Method 

For these systems we follow similar consideration as for the case of two-electron hydro- 
genic chains. The variational procedure is used to explore the problem. Physical relevance 
arguments are followed to choose the trial function (see, e.g. [10]) which is given by the 
function (7). 

n = 1. (aee) 

This case is only mentioned for the sake of completeness. It is known that the helium 
atom exists for any magnetic field strength [14]. At zero field and as well as for small 
magnetic fields B < 0.75 a.u. the spin-singlet state 1^0+ is the ground state. For 
B > 0.75 a.u. the spin-triplet state 1^(— l)""" becomes the ground state. Neutral 
Hehum atom is the least bound system among two-electron Coulomb systems made 
from a-particles. 

n — 2. (i) {aaee) 
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B (a.u.) 
System 


102 10^ 10^ 


4.414 X 10^3 G 


10^ 10^ 


10^ 


He 


25.65 54.37 106.4 


126.0 


191.4 319.7 


494.3 


He^+ 


33.98 80.49 174.51 


212.14 


343.47 616.68 


1016.75 


Hef 


26.58 68.93 163.90 


202.60 


352.50 684.19 


1212.40 


He6+ 






272.07 576.85 


1089.89 


HeH+ 


28.36 64.24 133.49 


160.50 


253.22 440.24 


709.65 


(H - He - H)2+ 


- 142.40 


172.58 


279.39 509.99 


843.38 


(He - H - He)^"^ 


- 153.62 


190.22 


320.63 603.91 


1029.95 


(H - He - He - H)^+ 






275. 585.0 


979.1 


(He - H - H - He)'^^ 






223. 510.4 


885.2 



TABLE VH: Double ionization energies Ei in Ry for the ground state ^H„ of the two-electron 
helium and helium-hydrogenic linear systems (finite chains) in a strong magnetic field (the ground 
state for (HeH)+ is ^H). 

The Heg"*" -molecule was studied in details in [15] in a magnetic field 5 = — 
4.414 X 10^^ G. It was shown that the lowest total energy state depends on the 
magnetic field strength. Similarly to the case of ppee, it evolves from the spin- 
singlet ^Tig metastable state at < 5 < 0.85 a.u. to a repulsive spin-triplet 
state (unbound state) for 0.85 a.u. ^ B < 1100 a.u. and, finally, to a 
strongly bound spin-triplet state. Hence, there exists quite large domain of 
magnetic fields where the Hc2^-molcculc is unbound being represented by two 
atomic helium ions in the same electron spin state but situated at the infinite 
distance from each other. The optimal geometry of the He2'''-molecule (when 
exists) corresponds always to the elongation along a magnetic line forming a 
finite chain. It is assumed that the chain in the state is stable towards 
the deviation from linearity. This chain is stable (or metastable) for all studied 
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S(a.u.) 

System 


102 


103 


10^ 


4.414 X 10^^ G 


10^ 


10^ 


10^ 


He2+ (R) 


0.463 


0.212 


0.106 


0.0902 


0.060 


0.0353 


0.023 


Hef {R,R) 


0.67, - 


0.27, - 


0.122, - 


0.116, - 


0.063, - 


0.0358, - 


0.023, - 


Rel+ {Ri,R2,Ri) 


— 






— 


0.089,0.060, - 


0.047, 0.037, - 


0.030,0.023, - 


HeH+ 


0.440 


0.203 


0.104 


0.092 


0.0585 


0.0356 


0.0238 


(H-He-H)2+ {Ri,Ri) 






0.105, - 


0.092, - 


0.059, - 


0.035, - 


0.022, - 


(He-H-He)3+ {Ri,Ri) 






0.095, - 


0.081, - 


0.051, - 


0.030, - 


0.018, - 


(H - He - He - H)'^^ 

{Ri, R2, Ri) 
(He - H - H - He)'^+ 

i?2, Rl) 










0.07, 0.10, - 
0.08, 0.12,- 


0.047, 0.030, - 
0.041, 0.025, - 


0.027, 0.015, - 
0.025, 0.019, - 



TABLE VIII: Internuclear equilibrium distances (in a.u.) for the ground state of the two- 
electron helium linear systems (finite chains) in a strong magnetic field (the ground state for 
(HeH)"^ is the ^H state). All configurations (except for (HeH)"^) have center of symmetry and 
symmetric interproton distances are not displayed. 

magnetic fields. However, this chain has the total energy higher than the Hcg"*" 
chain (see below) for S > 3 x 10^ a.u. and thus less preferable energetically. 

(ii) {apee) 

It is the simplest 2e mixed hehum-hydrogen system. A detailed study of the low- 
lying electronic states ^E, ^E, ^n, of the HeH"*" ion was carried out in [16]. The 
ground state evolves from the spin-singlet ^E state for small magnetic fields B < 
0.5 a.u. to the spin-triplet (unbound or weakly bound) state for intermediate 
fields and to the spin-triplet strongly bound state for B > 15 a.u. When 
the HeH"*" molecular ion exists, it is stable with respect to a dissociation. In 
the domain B > 15 a.u. the optimal geometry is linear and parallel: the ion is 
elongated along a magnetic line. Hence, the chain is formed. With a magnetic 
field increase the chain gets more bound and more compact. At magnetic fields 

24 



B < 10^ a.u. the double ionization energy Ej of the HeH"*" ion is smaller but 
comparable with one of the Hel'^ ion. However, for B > 10^ a.u. Ej gets, in fact, 
the smallest value among 2e Helium-contained molecular ions. 

(i) (aaaee) 

In field-free case the system (aaaee) does not display any binding. However, for 
magnetic fields B > 100 a.u. the He3'''-molecule becomes bound in the linear 
configuration aligned along the magnetic line. For 100 a.u. ^ B < 1000 a.u. the 
state is the ground state [18]. This state is a metastable state for any mag- 
netic field, its total energy lies above the total energies of its lowest dissociation 
channel. For B > 1000 a.u. the state ^H^ is the ground state. For magnetic 
fields 1000 a.u. ^ B <3 x 10^ a.u. the total energy of the dominant dissociation 
channel He^+ ^ Rel+i^U^) + a is lower than the total energy of the Hcg ion 
in the state. Thus, in this range of magnetic fields, the ion He3'''(^n„) is a 
metastable state towards the lowest channel of decay. Hence, for magnetic fields 
B >3x 10^ a.u. the molecular ion Heg"*" ion in the ^H„-state is stable. With an in- 
crease of the magnetic field strength, the total energy at the equilibrium position 
decreases, the system becomes more bound (in this case, the double ionization 
energy is Ei = —Et, it increases with B) and more compact (the internuclear 
equilibrium distance decreases with B). 

(ii) (papee) 

In field-free case the system (papee) does not display any binding. However, 
for magnetic fields B > 10^ a.u. the (H — He — H)^'^-ion becomes bound in the 
linear configuration aligned along the magnetic line with the ^H^ state as the 
ground state. This ion is stable. 

(iii) (apaee) 

In field-free case the system (apaee) docs not display any binding. However, for 
magnetic fields B > 10^ a.u. the (He — H — He)^^-ion becomes bound in the 
hnear configuration aligned along the magnetic line with the ^H„ state as the 
ground state. This ion is unstable towards a decay to He2''"(^H„) + p, however, 
at B > 10^ a.u. the ion (He — H — He)^"^ becomes stable. 
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(i) (4Q;2e) 

In field-free case the system (4Q;2e) does not display any binding. However, for 
magnetic fields B > 10^ a.u. the He^^-molecule becomes bound in the linear 
configuration aligned along the magnetic line with the state as the ground 
state. With an increase of the magnetic field strength, the total energy at the 
equilibrium position decreases, the system becomes more bound (in this case, the 
double ionization energy is Ej — —Et, it increases with B) and more compact 
(the internuclear equilibrium distance decreases with B). For magnetic fields 
B > 10^ a.u. the total energy of the dominant dissociation channel He3^(^n„)-|-Q; 
is lower than the total energy of the He4^(^n„) ion. Thus, the ion He4^('^n„) is a 
metastable state toward the lowest channel of decay. It is also unstable towards 
decay to He^+(^n„) + 2a for magnetic fields 10^ < S < 2 x 10^ a.u. 

(ii) {paap2e) 

In field-free case the system {paap2e) docs not display any binding. However, 
for magnetic fields B > 10^ a.u. the (H — He — He — H) ^'''-molecule becomes 
bound in the linear configuration aligned along the magnetic line with the ^H^ 
state as the ground state. With an increase of the magnetic field strength, the 
system becomes more bound (the double ionization energy increases with B) 
and more compact, i.e. both, the internuclear equilibrium distance Ri between 
a proton and the closest a particle, and the distance R2 between the two a 
particles, decrease with B. For magnetic fields B > 10^ a.u. the total energy 
of the dominant dissociation channel He2'''(^Hu) -|- 2p is lower than the total 
energy of the (H — He — He — H)^'*' ion. Thus, the ion (H — He — He — H)^"^ is a 
metastable state toward the lowest channel of decay. 

(iii) {appa2e) 

In field-free case the system {appa2e) does not display any binding. However, for 
magnetic fields B > 10^ a.u. the (He — H — H — He)^"''-molecular ion becomes 
bound in the linear configuration aligned along the magnetic line with the ^H„ 

state as the ground state. With an increase of the magnetic field strength, the 
system becomes more bound (the double ionization energy increases with B) and 
more compact, i.e. both, the internuclear equilibrium distance Ri between a 



26 



proton and the closest a particle, and the distance R2 between the two protons, 
decrease with B. For magnetic fields B > 10^ a.u. the total energy of the 
dominant dissociation channel He2^(^n„) + 2p is lower than the total energy of 
the (He — H ^ H — He) "'^-molecular ion, thus being a metastablc state toward 
the lowest channel of decay 

= 5. The results of the analysis of 5-center, 2-electron systems is shown in Table IX. It is 
not found an indication to binding of the proton-free systems (5a2e) for the whole 
domain of studied magnetic fields, while (4ap2e) gets bound at 5 ~ 10^ a.u. be- 
ing unstable decaying towards many different finite chains. The system {3a2p2e) is 
unbound although a particular configuration {apapa2e) displays a minimum in the 
potential curve. The two a-contained system are bound in both symmetric configura- 
tion - {papap2e) and {apppa2e) - while the latter one is more bound even for magnetic 
field B ~ 10^ a.u. This system is unstable with dominant decay mode to {apa2e). One 
Q!-contained system {ppapp2e) is bound at ~ lO'^a.u. and it is stable(!). It is worth 
noting that the system 5p2e is bound for magnetic fields B > 10^ a.u. (see Table III 
and a discussion on p. 15). 

= 6. It is not found an indication to the bound state of any 6-center system even for the 
highest magnetic field studied. 

C. Results 

The results of the calculations are presented in Tables VII-IX. Three traditional systems 
He, Hel'^ and HeH"^ exist for all studied magnetic fields 10^ < B < 10^ a.u. The first exotic 
molecular system Hef^ appears at ~ 100 a.u. in linear configuration and exists for larger 
magnetic fields. For 100 ^ B <5x 10^ a.u. the Hcg"^ ground state is a metastable state with 
respect to its lowest dissociation channel. For magnetic fields B > 5 x 10"' a.u. the ground 
state of the system Hep' becomes a strongly bound state. Another exotic molecular system 
He4''' appears at ~ 10^ a.u. as a metastable state. No other two-electron molecular helium 
systems are seen for B < 10^ a.u. At large magnetic field the ground state of all studied 
systems is the spin-triplet state with spin projection rus — —1 and total magnetic quantum 
number m — —1. For n > 1 the optimal geometry of the molecular system is linear, and the 
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Composition 


Configuration 


Et = -Ej (Ry) 
B = 10^ a.u. B = 10^ a.u. 


1-a 4-p 


H-H-He-H-H 


Unbound Bound 

~ -450 -866.0 

(Ri = 0.0228, = 0.0203 a.u.) 


2-a 3-p 


He-H-H-H-He 
H-He-H-He-H 


Bound Bound 

-414.5 -792.6 

Bound Bound 

-485.3 -873.9 

{Ri = 0.0306, i?2 = 0.0189 a.u.) 


3-a 2-p 


He-H-He-H-He 
H-He-He-He-H 


Unbound "Bound" 

~ -420 -860.0 

{Ri = 0.023, i?2 = 0.018 a.u.) 
Unbound Unbound 

~ -620 ~ -1055 


4-a 1-p 


He-He-H-He-He 


Unbound Bound 
~ -380 -862.4 

= 0.0356, R2 = 0.0195 a.u.) 



TABLE IX: 2-electron 5-center molecular ions (finite chains) in a magnetic field in state - 
symmetric, spin-triplet configuration parallel to the magnetic field direction. Et,Ej is total and 
double ionization energy, respectively. For unbound states a characteristic total energy indicated. 

28 



system is aligned along magnetic field. Thus, each molecular system forms a finite chain. It 
is checked that such a linear configuration is stable with respect to small vibrations and its 
vibrational energies can be calculated. However, we were not able to check stability of the 
configuration with respect to small deviations from linearity and to calculate the rotational 
energies. All studied finite chains are characterized by two features: with a magnetic field 
growth (i) their binding energies increase and (ii) their longitudinal lengths decrease - each 
system becomes more bound and compact. 

It is worth noting that among two-electron helium finite chains the system He^'^ in triplet 
^Uu state has the lowest total energy in the domain 10^ ^ B <3x 10^ a.u., whereas at larger 
magnetic fields 3 x 10^ ^ B ^ 10^ a.u. the finite chain Hcg^ in triplet "^Ilu state acquires 
the lowest total energy. In the domain 2 x 10^ ^ B < 10^ a.u. all studied two electron finite 
helium chains become stable with the only exception of He4'''. 

Conclusions 

A complete non-relativistic classification of one- two electron finite molecular chains (poly- 
mers) made out of protons/ a— particles in a strong magnetic field is presented. It is naturally 
assumed that the ground state of any one-electron chain is lag (or la for non-symmetric 
systems), while for any two-electron chain is spin-triplet ^11^ (or ^11 for non-symmetric sys- 
tems). All calculations were carried out in variational methods with state-of-the-art trial 
functions. Protons and a— particles are assumed to be infinitely-massive and situated along 
a magnetic line. 

It is clearly seen the existence of three magnetic field thresholds [22] , 

Bi^^ - lO^a.M. , ^ - W^a.u. , ^ - W^a.u. . 

At magnetic fields B < 10^ a.u. the only traditional ions, atoms and molecules may exist, the 
chains are not well-pronounced, they are very short containing at most two heavy particles. 
However, at 10^ < B < 10^ a.u. several new exotic ions appear in addition to traditional 
ones. All ions immediately form strongly bound linear chains aligned along a magnetic field. 
At S ~ 10^ a.u. several more new exotic ions appear quickly forming linear chains. Then 
similar appearance of new exotic ions happens at S ~ 10^ a.u. It is quite interesting that 
the ions which already appeared (existed) below some magnetic field threshold, above of the 
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threshold they become stable. It is worth noting that for fixed magnetic field the neutral 
systems are always the least bound ones. 

Concluding we present a list of l-2e proton- aparticle contained ions for which a certain 
magnetic fields exist where they are stable, 

H, H+, H2+,He+, He3+, (HeH)^^ , 

H- , H2 , H+ , , E.I+ , He , He+ , He^+ , (HeH)+ , (HHeH)^+ , (HeHHe)^+ , (HHHeHH)^" 

among the 25 Coulomb l-2e systems which (can) exist in a magnetic field (see Tables I- IX) . 

All presented results are obtained in non-relativistic way with an assumption that masses 
of heavy particles are infinite. They can be considered as an indication to a new atom- 
molecular physics in magnetic fields B > 10^ a.u. It encourages us to an exploration of finite 
mass effects in a magnetic field. This issue looks quite complicated technically due to absence 
of a separation of variables, especially, in the case of more than two particles and non-zero 
total charge of the system. Those two cases are exactly ones which are the most important 
from the point of view of obtained results: the most bound systems contain usually more 
than two bodies and charged. Another important issue is related to relativistic corrections to 
our non-relativistic results. Although in our understanding the Duncan qualitative argument 
[9] sounds physically, it needs to be checked quantitatively. Present authors plan to study 
both issues in near future. 
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